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Using techniques from symplectic geometry, we determine when a pure state of three qubits is 
up to local unitaries uniquely determined by its reduced density matrices. We moreover show 
that this is always the case if one is given the additional promise that the quantum state is not 
convertible to the Greenberger-Horne-Zeilinger (GHZ) state by stochastic local operations and 
classical communication (SLOCC). 



1. Introduction 

Finding solutions to Hamilton's equations for a given system is a standard problem in 
classical mechanics. The system in question is often invariant with respect to a certain 
group of symmetries, and when the symmetries are continuous, such as in the case of 
rotational symmetry, they form a Lie group K. It is well understood nowadays that 
the existence of symmetries is inevitably connected to the first integrals of Hamilton's 
equations. The celebrated theorem of Arnold (1989) states that in case when K = 
is an n-dimensional torus and n half the dimension of the phase space Ai0 then the 
system is completely integrable. When the group K is not abelian, it is still possible 
to find corresponding first integrals, although they typically do not Poisson commute. 
For example, when a particle is moving in a potential with rotational symmetry, so 
that K = SO (3), then the conserved quantities are the three components of the angular 
momentum, corresponding to the invariance of the system with respect to infinitesimal 
rotations about the axes x, y, and z. These infinitesimal rotations generate the Lie algebra 
SO (3) of 50(3). There are of course many possible bases of so (3), each corresponding to 
a different choice of rotation axes and each giving three first integrals. The mathematical 
object which encodes information about the first integrals for all possible choices of 
generators of so(3) is an equivariant map fi : M— )-so(3)* from the phase space M to the 
space of linear functionals on the Lie algebra so (3). For every infinitesimal symmetry 
^€so(3) one obtains a corresponding first integral by the formula fi^{x) = {fi{x) , ^), 
where by (•,•) we denote the pairing between linear functionals from 6* and vectors 
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^ A phase space is a symplectic manifold {M,ui), where is a closed nondegenerate two-form. 
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in t. This idea can be generalized to arbitrary Lie groups K and a corresponding map 
;U : M — )• is called a momentum map (Guillemin &; Sternberg 1990). 

Remarkably, momentum maps appear naturally not only in classical but also in 
quantum mechanics. Indeed, the Hilbert space Ti on which a given quantum-mechanical 
system is modeled can be seen as a phase space if we identify vectors that differ by a global 
rescaling or a phase factor e^"^. The set of (pure) quantum states is thus isomorphic to the 
complex projective space P('H), which is well-known to be a symplectic manifold. When 
the considered system consists of subsystem^ then the space Ti has the additional 
structure of a tensor product, namely Ti =71^^, where Tii is the single-particle Hilbert 
space. The mathematical properties of this structure manifest physically as entanglement, 
i.e. quantum correlations between subsystems. These correlations are invariant with 
respect to the local unitary action, i.e. to the action of the Lie group K = SU{'Hi)^^ on 
Ti by the tensor product. Since K preserves the phase space structure of P('H) one gets a 
momentum map // : P(T-L) — )• , and the image //([f]) of a quantum state [v] G P(^) is, up 
to some unimportant shifting, the collection of its one-body reduced density matrices (see 
Sectional). The matrices pi can be diagonalized by the action of K and their eigenvalues 
can be ordered, for example decreasingly. The map which assigns to the state [v] the 
collection of the ordered spectra of its one-body reduced density matrices will be denoted 
by ^ : P(^) — )■ t^. This definition can be generalized to an arbitrary compact Lie group 
K, with the positive Weyl chamber corresponding to a choice of a Borel subgroup 
B cG of the complexification G = K'^. 

In the early '80s, the convexity properties of the image of momentum map, n{M), 
were investigated by Atiyah (1982) and Guillemin and Sternberg (1982, 1984) in the case 
of the abelian group K = T". They proved that fJ,{M) is a convex polytope; it is in fact the 
convex hull of the image of the set of fixed points. For nonabelian K, this is typically no 
longer true. However, as it was shown by Guillemin &: Sternberg for Kahler manifolds and 
by Kirwan (1982, 1984) for general symplectic manifolds, the image ^(M) = p-{M) R 
is always a convex polytope — the so called Kirwan polytope^ or momentum polytope. 
Finding an explicit description of this polytope is in general a difficult problem; it is an 
instance of the quantum marginal problem, or A^-representability problem in the case of 
fermions (Ruskai 1969, Coleman & Yukalov 2000). Following the groundbreaking work of 
Klyachko (1998) on the famous Horn's problem, which can be formulated as the problem 
of determining a Kirwan polytope, Berenstein & Sjamaar (2000) found a general solution 
for the case where M is a coadjoint orbit of a larger group. This was in turn used to 
solve the one-body quantum marginal problem; namely, sets of inequalities describing the 
Kirwan polytope ^(P('H)) were given (Klyachko 2004, Daftuar & Hayden 2004, Klyachko 
2006). Intriguingly, the Kirwan polytope can under certain assumptions be also described 
in terms of coinvariants or, equivalently, in terms of the representations that occur in 
the ring of polynomial functions on M (Ness 1984, Brion 1987). In the context of the 
quantum marginal problem, this has been discovered by Christandl et al. (2006, 2007) 
and Klyachko (2004) . The interplay of these two complementary perspectives can also be 
seen in (Christandl et al. 2012), where an algorithm has been given for the more general 
problem of computing the distribution of eigenvalues of the reduced density matries of 
a multipartite pure state in P('H) drawn at random according to the Haar measure; in 
particular, this gives an alternative solution to the one-body quantum marginal problem. 

The above line of work can also be seen as one of the first successful applications 
of momentum map geometry to the theory of entanglement (see also Klyachko 2008). 



^ We assume for simplicity that their Hilbert spaces are of the same dimension. 
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In 2011, the importance of the momentum map to entanglement was investigated from 
a different perspective by Sawicki et al. The authors showed that restricting the map 
^ to different local unitary orbits in P(?^) gives rise to a well defined purely geometric 
measure of entanglement. They also pointed out that the properties of the fibers of ^ are 
crucial to the solution of the local unitary equivalence problem and gave an algorithm for 
checking it for three qubits (Sawicki & Kus 2011). Geometrically, the problem of local 
unitary equivalence of two states reduces to determination whether they both belong to 
the same orbit of the local unitary group K = SU {Hi}^^ . Form the physical point of 
view it corresponds to checking if one of the states can be obtained from the second one 
by non-dissipative quantum operations restricted to subsystems - an important problem 
in quantum engineering of states aiming at practical applications (Kraus 2010), (Kraus 
2010a). 

In the subsequent paper the authors analyzed the geometric structure of the fibers 
of for two distinguishable particles, two fermions and two bosons (Huckleberry et 
al. 2012). In all these cases, the JT-action is spherical, i.e. the Borel subgroup B C K'^ 
has an open orbit in By Brion's (1987) theorem, this implies that every fiber of ^ 

is contained in a single J^-orbit. That is, in these cases every quantum state is up to local 
unitaries uniquely determined by the spectra of its reduced density matrices. Moreover, 
each fiber of ^ has the structure of a symmetric space (Huckleberry et al. 2012). 

In situations involving larger numbers of particles, e.g., iV-qubit systems with N >2, 
the action of K is not spherical and Brion's theorem cannot be applied. The identification 
of i^-orbits, i.e. classes of states which can be mutually transformed into each other by 
local unitary transformations, generically requires more information than it is contained 
in the spectra of the single-particle reduced matrices. The main goal of this paper is to 
investigate the set of quantum states which are mapped by ^ to the same point of the 
Kirwan polytope ^(M) (again, these are quantum states whose reduced density matrices 
have the same spectra) by studying the fibers of ^ or, equivalently, the symplectic 
quotients Ma = '^~^{a)/K, in the case where the action of K is no longer spherical. 
Specifically, we consider the above-mentioned case of qubits, where the Kirwan 
polytope is known explicitly (Higuchi et al. 2003, Bravyi 2004). A detailed analysis is 
carried out for three qubits and we make some remarks about the general case. Our main 
tools are the convexity theorem for projective subvarieties (Ness 1984, Brion 1987), as 
well as some general properties of orbit spaces which we recall in Section [2j Specifically, 
we show that Mq, is generically a two-dimensional stratified symplectic space (Sjamaar & 
Lerman 1991). For the points in the boundary of the Kirwan polytope the situation is 
very different. We prove that in this case Mq, is a single point, i.e. the dimension drops 
down by two compared with the interior. In particular, these results characterize the 
ii'-orbits which are uniquely determined by the spectra of the reduced density matrices. 
They may be contrasted with the well-known fact that the two-particle reduced density 
matrices generically suffice to determine a pure-state of three qubits (Linden et al. 2002). 

The complexified group G = K'^ plays its own role in the classification of states of 
multiparticle quantum systems. Its elements correspond to stochastic local operations 
with classical communication (SLOCC); see, e.g. Horodecki et al. (2009). States can be 
classified by identifying those which belong to the same SLOCC class, i.e. the same G- 
orbit. As in the case of locally unitary equivalent states (orbits of K), a detailed analysis 
of the corresponding Kirwan polytopes is useful in such a classification. In this paper we 
therefore examine the Kirwan polytopes 'l/(G.[f]) for the three-qubit SLOCC classes (Dur 
et al. 2000), i.e. for the orbit closures of the complexified group G = K'^, and describe 
their mutual relations. In particular, we find that the map ^ separates iC-orbits when 
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restricted to the closure of the so-called Vl^-class. That is, states from the Vl^-class are 
up to local unitaries characterized by the collection of spectra of their one-qubit reduced 
density matrices. 

The paper is organized as follows. In Section [2] we gather useful facts about orbit 
spaces, especially concerning the principal and maximally-dimensional orbits. In Section[3] 
we recall the notion of a momentum map and state precisely various versions of the 
convexity theorems. In Section |4] we discuss in detail the structure of the fibers of ^ 
for three qubits. Throughout the paper we prove only new results and otherwise give 
references to the literature. 



2. Qualitative properties of orbit spaces 

In this section we gather some useful facts about orbits of group action on a connected 
manifold M. A good reference is (Bredon 1972). 

Let K he a, compact connected Lie group acting smoothly on a connected manifold 
M (dim Af = n), i.e. we have a map ^ : K x M — > M, ^{g, x) = ^g{x), such that 

^ghix) = <^gi<l>hix)) VxEM, (2.1) 

^eix) = x VxGM. 

Any orbit K.x = ^{K,x) of the i^-action is isomorphic to ^/Kx, where Kx C -ftT is the 
isotropy subgroup of x, i.e. the closed subgroup consisting of the group elements fixing 
X. We say that two orbits are of the same type, 

type(^M) = type(^M), (2.2) 

if and only if Hi and H2 are conjugate in K. It is easy to see that (|2.2p defines an 
equivalence relation and hence divides the space of i^'-orbits into classes. It is also possible 
to introduce a partial order on this set, and hence on the set of orbit types, by imposing 
that 

type{K/H,)>type{K/H2) 

if and only if Hi is conjugate to a subgroup of H2. Slightly abusing notation, we will 
write k/h for the corresponding orbit type. The following theorem is an important result 
for the classification of orbit types (Bredon 1972): 

Fact 1 (Principal orbit type). There exists a single maximal orbit type ^/h for the 
K -action on M , i.e. a subgroup H which is conjugate to subgroups of all other isotropy 
subgroups. Moreover, the union M(^f{-j of the orbits of type ^/h is connected, open and 
dense in M . 

The subgroup H C K described by Theorem [1] is called the principal isotropy subgroup, 
and the orbits in Mf^j^-^ are called principal orbits. Other orbits are classified as follows: 

• an orbit of type ^/Hi is called singular if and only if dim(^^i/_f/) > 0. 

• an orbit of type ^/Hi is called exceptional if and only if dim(^i/_ff) = and Hi^h is 
a non-trivial finite group. 

Let us for the remainder of this article assume that there exists a point x £ M such that 
the isotropy subgroup is discrete. By Theorem [H the principal isotropy subgroup is 
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conjugate to a subgroup of K^, hence it is also discrete. So if one point xGM has a 
discrete isotropy subgroup then there is an open dense connected set of points with a 
discrete isotropy subgroup. In fact (Montgomery et al. 1956): 

Fact 2. The union Mmax of the orbits of maximal dimension, i.e. the union of the 
principal and exceptional orbits, is a connected open dense set whose complement has 
dimension at most n — 2. If in addition M is compact, then there are only finitely many 
orbit types (i.e. finitely many different isotropy subgroups up to conjugacy). 

The following theorem gives an estimate for the dimension of the sets of singular and 
exceptional orbits (Montgomery et al. 1956): 

Fact 3. Let r = minjdim Kx : x G M}, i. e. r is the dimension of the principal isotropy 
subgroup. Let Mi = {x £ M : dimKx > r + i} where i>l. Then Mi is a subset of M 
of codimension at least i + 1. If M is compact and Hi{M, Z2) =0 then singular and 
exceptional orbits form a closed set of dimension at most n — 2. 

In sections [3] and |4] we will see that the existence of the principal orbit type is crucial to 
the classification of the fibers of the momentum map. 

3. Momentum map 

Let be a compact connected Lie group acting on a symplectic manifold (M, uj) in a 
symplectic way, i.e. the map $g from (|2.ip satisfies the additional condition that <I>*aj = u 
for any g £ K, where <I>* stands for the puUback of the form by ^g. Denote by T{M) the 
space of smooth functions on M. The symplectic structure on M induces the structure 
of a Lie algebra on J^(M), with the Lie bracket given by the Poisson bracket, 

{/, g]=u:{Xf,Xg), 

where Xf, Xg are the Hamiltonian vector fields associated to f,g£j^{M), i.e. df = 
uj{Xf,-) and dg = uj{Xg, ■). Hamiltonian actions are defined in such a way that every 
fundamental vector field arises from a Hamilton function in a consistent way (Guillemin & 
Sternberg 1990): 

Definition 1. We say that a symplectic action of K on (M, w) is Hamiltonian if and 
only if there exists a momentum map fi: M ^t* , i.e. a map which satisfies the following 
three conditions: 

1. For any ^ G i, the fundamental vector field ^{x) = J^|j_g$(e*^, x) is the Hamiltonian 
vector field for the Hamilton function fi^{x) = {fi{x) , i.e. dfj,^=u}{£,,-). 

2. The induced map t 9 ^ 1— t- /x^ G ^ homomorphism of Lie algebras, i. e. 

3. The map is equivariant, i.e. /^(^^(x)) = Ad* /i(x), where Ad* is the coadjoint 
action of K on t* defined by (Ad* a, S,) = (a, Ad^-i ^) in terms of the adjoint action 
Adg^= ^1^=05 6*^ of K on its Lie algebra t. 

For semisimple K, hence in particular for K = SlKj-Li)^^ , the momentum map fi is 
uniquely defined by the above properties (Guillemin & Sternberg 1990). 
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(a) Convexity properties of the momentum map 

We will now assume that M is compact and connected. Let us choose a maximal 
torus T C K, with Lie algebra t, and a positive Weyl chamber C t*. Denote by : M — )• 
fi{M) n the map which assigns to x £ M the unique point of intersection n{K.x) t^. 
Then the following convexity results hold: 

1. The image ^'(M) is a convex polytope, the so-called Kirwan polytope (Guillemin & 
Sternberg (1982, 1984) and Kirwan (1982, 1984)). 

2. The fibers of ^ (and hence the fibers of ^) are connected (Kirwan (1982, 1984)). 

3. The map ^ is an open map onto its image, i.e. for any open subset U C M the 
image "^{U) is open in ^'(M) = ;u(M) n t+ (Knop 2002). 

4. The image ^'(Afjnax) is convex (Heinzner &; Huckleberry 1996). 

Denote by G = K"^ the complexification of K. This is a complex reductive group; in 
particular, a choice of positive Weyl chamber is equivalent to a choice of Borel 
subgroup B C G. Let us assume as in (Brion 1987) that M C P('H) is a non-singular G- 
invariant irreducible subvariety of the complex projective space associated with a rational 
G-representation Ti, and let us choose a X-invariant inner product on Ti. Then M is a 
symplectic manifold when equipped with the restriction of the Fubini-Study form, 

a;[.](iH, A.]) = 2Imi41^ = -i^^^j^ VA, B G u(^), (3.1) 

where [v] G ¥(T-L) is the projection of a vector v £71, and where A is the fundamental 
vector field generated by the action of A G u{l-L). More concretely, we can represent the 
tangent space T[^,]P(?^) by T-L/Cv = (Cv)-*". In this picture, the tangent vector A^, is given 

by the orthogonal projection of Av onto (Cv)"*". M also carries a canonical momentum 
map cP('H)^r 

(MM),A) = -iMM vAEt, (3.2) 



In this situation, there are further convexity results for the closures G.x and B.x of G- 
and i?-orbits, respectively (which in general are singular varieties for which the above 
results do not apply): 



(v) The image ^{G.x) is a convex polytope (Brion 1987, Ness 1984). 

(vi) The image ^{B.x) is a convex polytope (Guillemin & Sjamaar 2006). 

(vii) There is an open dense set U cG.x such that ^(B.y) = ^{G.x) for y£U 
(Guillemin & Sjamaar 2006). 

(viii) There is an open dense set U C M such that ^'(M) = ^{B.x) = ^{G.x) for x £U 
(Guillemin & Sjamaar 2006). 

(ix) The collection of different polytopes ^(B.x), where x ranges over M, is finite 
(Guillemin & Sjamaar 2006). 
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Since Mmax is a connected, open and dense subset of M, (iv) implies the following: 

Fact 4. The set /i(Mjnax) H is an open, dense, connected, convex subset of^{M) = 
IJi{M) n t^. In particular, it contains the (relative) interior of the Kirwan polytope. 

(b) Fibers of the momentum map 

The following two facts characterizing kernel and image of the differential of the 
momentum map are well-known consequences of the definition (Guillemin &; Sternberg 
1982): 

Fact 5. The kernel of d^\^:TxM ^t* is equal to the (jJ -orthogonal complement of 
Tx{K.x), i.e. 

KeT{d^l\J = {Ye T^M : u{i F) = G = (T.(i^.x))^" . 

Fact 6. The image ofdfi\^ : TxM — t- i* is equal to the annihilator oft^, the Lie algebra 
of the isotropy subgroup C K . 

These results have immediate consequences for the characterization of the fibers of //: 
Using Fact [6l we observe that d^ is surjective if and only if the Lie algebra tx is trivial. 
This in turn happens if and only if the isotropy subgroup Kx is discrete. Therefore, if 
d/i is surjective at a single point x £ M (as we have assumed above) then by Fact [2] it is 
automatically so for all points in Mmax = {x £ M : dimii'.rE is maximal}. In particular, 
the implicit function theorem together with Fact [5] implies that for all x £ Mmax) the 
tangent space at x of the ^u-fiber through x has dimension 

dimTxifi-^ifiix))) = dim {Tx{K.x))^'^ = diniM - dim!* (3.3) 

and hence: 

Fact 7. For the points xSMmax; for an open, connected and dense subset 
of M, the intersection /i~"^(//(a;)) Pi Mmax of the fi-fiber through x with Mmax is a 
dim{Tx{K.x)) -dimensional manifold. 

Remark 1. An analogue of L'J. 'J\) is true in the case when the principal isotropy 
subgroup is not discrete, as can be confirmed by using the constant rank theorem instead 
of implicit function theorem. 

Understanding the fibers /i~^(a) away from the regular points in Mmax is a more delicate 
problem. Since the fibers of /x are X-invariant, it is convenient to introduce the symplectic 
quotient Mq, := (a) / Ka = (a) / K , which is in general a stratified symplectic space 
(Sjamaar & Lerman 1991). It is moreover connected by property (iii) in Subsection [a) 
The following is then immediate: 

Fact 8. The fiber f^i~^{a) intersects a single K -orbit if and only if Ma is a single 
point or, equivalently, if and only if is zero-dimensional. 

If a G ^'(Minax) then the maximal-dimensional stratum of Ma is simply {fj,~^{a) Pi 
MraoK) / Ka (Meinrenken &: Woodward 1999). Since the isotropy subgroup of any point 
in Mjnax is discrete, 

dim Ma = dim M - dim T - dim Ka (3.4) 

In other words, we can by mere dimension counting determine whether a X-orbit in 
Afmax is uniquely determined by its image under the momentum map. 
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4. Fibers of the momentum map for three qubits 

With the momentum map machinery presented in the preceding section we are now weh- 
equipped to analyze when a pure state of three qubits is, up to local unitaries, determined 
by the spectra of its reduced density matrices. 

To this end, let M = ¥{T-L) be the projective space of pure states associated with the 
three-qubit Hilbert space = (g) (g) C^. The group of local unitaries K = SU{2)^^ 
and its complexification G = SL{2)^^ act on M by the tensor products. The Lie algebra 
of K is t = 5u(2) © su(2) © su(2), i.e. triples of traceless antihermitian matrices. As 
described in the preceding section, M is a symplectic manifold with respect to the 
Fubini-Study form (|3.ip and the i^-action is Hamiltonian with canonical momentum 
map (ITD . 

It is easy to see that under the identification of t* with t induced by the trace inner 
product, the image /^([f]) is given by the collection of one-qubit reduced density matrices, 
namely 

KM)=iiPi-ll,P2-ll,P3-ll), (4.1) 



where / is the 2x2 identity matrix (see, e.g., (Sawicki et al. 2011, Christandl et al. 2012)). 

Let us fix the maximal torus T C to be the set of unitary diagonal matrices with 
determinant equal to one. Then the Lie algebra t is equal to the space of traceless 
antihermitian diagonal matrices. We choose as the positive Weyl chamber the following 
set of matrices: 

t+ = {(diag(-'jAi,iAi), diag(-iAi, iAi), diag(-iAi, zAi)) : Aj > 0} . (4.2) 

It follows that, up to some rescaling and shifting, the map ^ sends a pure state [v] to 
the (ordered) spectra of its reduced density matrices pi,p2,P3- 

The following theorem describes the Kirwan polytope in terms of inequalities for a 
general system of qubits: 



Fact 9 ((Higuchi et al. 2003)). For a N-qubit system, the constraints on the one-qubit 
reduced density matrices pi of a pure state are given by the polygonal inequalities 

where Pi^\ denotes the minimal eigenvalue of pi (i = 1, . . . , N ). 



The Kirwan polytope for three qubits is shown in Figure [TJ 
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Figure 1. The Kirwan polytope '^/{M) = ii{M) n t+ for three qubits. 



This polytope has 5 vertices, 9 edges and 6 faces. As a convex set it is of course generated 
by its vertices, which are 

t^SEP = {diag(0, 1), diag(0, 1), diag(0, 1)}, 
VBl = {diag(0, 1), diag(V2, 1/2), diag(i/2, 1/2)}, 
VB2 = {diag(V2, 1/2), diag(0, 1), , diag(V2, 1/2)}, (4.3) 
VB3 = {diag(V2, 1/2), diag(V2, 1/2), diag(0, 1)}, 
VGHZ = {diag(V2, 1/2), diag(V2, 1/2), diag(V2, 1/2)}- 

Since the polytope is full-dimensional, Sard's theorem (Sternberg 1964) implies the 
existence of a regular point in M, i.e. a point with discrete isotropy subgroup. Hence: 

Lemma 1. The set Mmax C M is connected, open and dense and consists of orbits of 
dimension dim K = 9. 

We will now analyze points inside the interior of the Kirwan polytope. Notice first that 
by Fact [3]the the preimage of any such point a contains a point x € Mmax- Therefore, 
Fact [7] and ()3.3p show that /i~^(a) n Mmax is a manifold of dimension 

dim(/i~^(a) n M^ax) = dim(r^K.x)-^'^ = dimM - dimt* = 14 - 9 = 5. 

Since Ka = T for points in the interior of the positive Weyl chamber, this manifold 
consists of 3-dimensional T-orbits, and (j3.4p implies that 

dimMa = 5-3 = 2. 

If we replace n~^{x) by ^'^^(x) = K.fi^^{x), we replace each T-orbit by a i('-orbit and 
therefore increase the dimension by 

dim Qq, = dim ^/Ka = dim K — dim T = 6, 

the dimension of the corresponding coadjoint orbit Qa = K.a C t* . The general formula 
for dimension of dimOo, is 

Ki K2 Ka 

dim 17„ = dim K - (^ „ + ^ + ^ - 3) , 

n=0 n=0 n=0 
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where Ki is the number of distinct eigenvalues of pi and TfLi 'fi the niultiphcity of Tz-th 
distinct eigenvalue of pi (Sawicki et al. 2011). Summing up, we proved: 

Theorem 1. For any point a inside the interior of the Kirwan polytope there exists 
a point x G Mmax such that a = p{x). The manifold H Mmax is b- dimensional, 

consisting of 3-dimensional T -orbits. Moreover, ^'"^(q;) H Mmax is an 11 -dimensional 
manifold consisting of 9 -dimensional orbits K.y with the property p{K.y) = ^a> and the 
symplectic quotient has dimension 2. 

The following is a direct consequence of Theorem [1] (cf. the discussion at the end of 
Section [3]) : 

Corollary 1. The orbits K.x for which p{K.x)r\^j^ belongs to the interior of 
Kirwan polytope cannot be separated by the momentum map. In other words, a pure state 
of three qubits whose spectrum is non-degenerate and satisfies the polygonal inequalities 
with strict inequality is never determined up to local unitaries by the spectra of its reduced 
density matrices. 

What is left is to analyze points in the boundary of the Kirwan polytope. We postpone 
this to the end of the section (see Subsection [cj) . 

(a) The SLOCC classes and their Kirwan polytopes 

As mentioned in the Introduction, the classification of G-orbits, where G = K^ 
gives another view on the entanglement properties of composite system states. In this 
subsection we explicitly compute the Kirwan polytopes ^{G.x) for all G-orbit closures, 
and show how they are related to the polytope ^{FiT-L)). Physically, G-orbits correspond 
to SLOCC classes (Dur et al. 2000), hence our results describe the spectra of the reduced 
density matrices of the quantum states in each SLOCC entanglement class. The problem 
of classifying G-orbits in is inherently connected to the momentum map geometry as 
well as to the construction of the so-called Mumford quotient (Kirwan 1984, Ness 1984). 
We explain this connection in (Sawicki et al. 2012, Walter et al. 2012). 

Classification of G-orbits It has been shown in (Dur et al. 2000) that there are six 
SLOCC entanglement classes, i.e. G-orbits in P('H). For convenience of the reader we list 
them below and briefly summarize their basic geometric properties: 

1. The G-orbit of the Greenberger-Horne-Zeilinger state, 

xqbz = ["^ (|000) + |111))] (Greenberger et al. 1989): It is the open dense orbit, 
hence of (real) dimension dimP('H) = 14. This also follows from the facts that 
/i(i^.XGHz) = and that the orbit K.xquz is Lagrangian, hence of dimension 7, as 
was shown by Sawicki & Kus (2011), since for any i^-orbit in /i~^(0) we have the 
relation dimG.x = 2 • dimi^.x. 

2. The G-orbit of the M^-state, xw = (1 100) + |010) + |001))]: Here, dimG.xvK = 
12, while dim K.xw =8. For the proof, it is enough to compute the dimension 
of the tangent space Txy^riG.xw), which can be represented as the projection of 
Span^lArEvy : Ag q} onto the orthogonal complement of xw- The Lie algebra q is 
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equal to 512(C) s[2(C) s[2(C), where 

s[2(C) = Spanc {E12, E2U Eu - ^22} , 

and Eij is the 2x2 matrix with a single non-zero entry equal to one in the i-th. 
row and j-th. column. It is easy to see that the following seven vectors 

{E12 / I)xw = (/ E12 I)xw = (/ / Ei2)xw oc |000), 

{E21 0/0 I)xw oc 1 110) + 1 101) , 

(/ E21 I)xw oc |110) + |011), 

(/0/0^2l)xH/OC |101) + |011), 

((^11 -£;22)0/0/)xH'Oc-|lOO) + |010) + |001), 

(/ (^11 - E22) I)xw oc +|100) - |010) + |001), 

(/ / (^11 - E22))xw oc +|100) + |010) - |001) 

span a complex vector space of dimension 6 after projection onto (the last 
3 vectors become linearly dependent). We conclude that dim G.xvi/ = 2 • 6 = 12. 
Similarly, one shows that dimK.xw = 8. 

3. The G-orbits through the bi-separable Beh states xbi = [^(|000) + |011))] = [|0) 
^(|00) + |11))], XB2 = [^(|000) +I101))], and = ([;^(|000) + |110))]: In a 
similar way as for xw one can show by an explicit computation that dim G.XBk = 8 
and dim. K.XBk = 5. 

4. The G-orbit of separable states, generated by xsep = [|000)]: By the Kostant- 
Sternberg theorem, it is the unique symplectic orbit (Sawicki et al. 2011). In fact, 
since |000) is in the ii'-orbit of the highest weight vector, the Borel-Weil theorem 
(Sepanski 1995) implies that G. xsep = K. x^ep — K/T = G/B; in particular, the 
orbit is Kahler and dimG.xsEP = dimK.xsEP = 9 — 3 = 6. 



Kirwan polytopes of G-orbit closures Let us write Xj := G ■ xj for the closures of 
these G-orbits. The representatives Xj that we have chosen above satisfy the following 
property: 

Lemma 2. We have '^{xj) = fJ.{xj) =Vj, where Vj is defined in lji4-3\ ). Moreover, ^{xj) 
is the closest point to the origin in the Kirwan poly tope ^(Xj). 

The first claim follows from a simple computation. The fact that vj is the closest point 
to the origin of the Kirwan polytope of the corresponding G-orbit closure follows from 
the general theory of (Kirwan 1982): the gradient descent with respect to the norm 
square of the moment map is at any point x implemented by the vector field generated 

by ifi{x) (zHQq, and one easily checks that i/x(xj)|_^ =0. It will be explained in more 

detail in (Sawicki et al. 2012, Waher et al. 2012). 

We will now describe the Kirwan polytopes of the G-orbit closures: Since G.xghz is 
dense in P(^), we immediately obtain that ^(Xghz) is equal to the full Kirwan polytope 
for three qubits as described by Fact [9l On the other hand, since G.xsep = -f^-^^SEP is 
a single X-orbit, ^'(XgEp) is a single point, namely ^'(^sEp) = {'WSEp}- Similarly, one 



12 



A Sawicki, M Walter, and M Kus 



finds that the Kirwan polytopes for the bi-separable Bell states XBk are equal to the 
one-dimensional line segments from VBk to vsep- Therefore, the only non-trivial task is 
the computation of the Kirwan polytope for the M^-state: 

Proposition 1. The Kirwan polytope ^(Xiv) is equal to the convex hull of the points 
vbi, vb2, vb3 and vsep- In particular, it is of maximal dimension. 

Proof. Clearly, ^(Xw) is a subset of the full Kirwan polytope ^{¥{T-L)). On the other 
hand. Lemma [2] and convexity imply that it is also contained in the half-space through 
vw with normal vector vw- Since the intersection of this half-space with ^(F('H)) is 
precisely equal to the convex hull of the points vbi, vb2, vb3 and fsEP) we only need to 
show that these points are contained in the Kirwan polytope. 

We will in fact show that the corresponding preimages Xj are contained in the orbit 
closure Xyy. For this, we observe that the action of the complexification C G of the 
maximal torus T C K applied to xw gives rise to all states of the form 

[ci|100)+C2|010) + C3|001)] (c,/0) (4.4) 
In particular, xj C Xw = G ■ xw for j = Bl, B2, B3, SEP, and the claim follows. ■ 

(b) Sphericality of the W SLOCC class 

In this subsection, we show that Xyy = G.xw is a spherical variety. It then follows by 
Brion's theorem that every quantum state in Xw is, up to local unitaries, characterized 
by the collection of the spectra of its one-qubit reduced density matrices. In other words, 
^ separates the ii'-orbits in X^/. 

We start by clarifying the geometric structure of X\y. Note that since G = is a 
complex reductive group, the closures of any orbit closure Xj is a G- invariant irreducible 
subvarieties of P(^); however, these varieties will in general be singular. This is in fact 
already the case for Xw. Indeed, it is known from (Klyachko 2008) that 

Xw = P(H) \ (G.XGHz) = {H G IP(^) : Det(^;) = 0}, 

where Det is the Cayley hyperdeterminant (the basic invariant for the G - representation 
%), and one readily verifies that the tangent space at xsep £ has complex dimension 
7 > 6 = dime Xw ■ 

Let us denote by ^w '■ ~^ ^* the restriction of : P('H) — t- t* to Xw- Our aim now 
is to prove that fiw separates the if-orbits in Xw- To this end we use Brion's theorem 
(cf. (Huckleberry & Wurzbacher 1990)): 

Fact 10. (Brion 1987) Let G = K'^ be a connected complex reductive group, 
% a rational G -representation, and X a G-invariant irreducible subvariety of V{T-L) 
(cf. Subsection\W^ . Then the following are equivalent: 

1. X is spherical, i.e. the (every) Borel subgroup B has a Zariski-open orbit in X. 

2. For every x£X the fiber /_f~^(/_f(x)) is contained in a single K-orbit, K.x. 
We will now show that indeed: 

Proposition 2. The G -variety Xw is spherical. 
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Proof. Consider the Borel subgroup B, which consists of the lower-triangular matrices 
in G. Its Lie algebra b is equal to © n, where 

n = Spanc{-E2i ® (g) /,/ ® / ® ^21}- 

In order to show that Xy\r is spherical, we have to show that there exists a Zariski-open 
orbit B.x. Since B.x is Zariski-open in its closure, it suffices to show that the closure 
of B.x is equal to Xy\r. Now, since the closure of B.x is a closed subvariety, it is either 
equal to Xw or of lower dimension. Therefore, it suffices to show that the dimension of 
B.x is equal to the dimension of Xy[/. Since B.x is a smooth variety, we can compute 
its dimension by computing the dimension of the tangent spaces at any point. Hence in 
order to show that Xw is spherical, it suffices to show that dime Tx{B.x) = dime X\\r = 6 
for any single point x in the G-orbit of xw- We will consider the state 

x = [-^(|100) + |010) + |001) + |000))] = ((^ J \ ^(^I®I)-xw. 

Indeed, one can easily verify that the tangent vectors generated by b, i.e. the projection 
of 

{E2i®I®I)xo^\im) + |101) + |100), 
(7 ® E21 18) I)x oc |110) + |011) + |010), 
(/(8)I(g)^2i)2;oc |101) + |011) + |001), 
((^11 - E22) (^I^I)x oc -|100) + |010) + |001) + |000), 
(/ ® {En - E22) I)x oc +|100) - |010) + |001) + |000), 
{I®I®{Eii-E22))xo^+\im) + |010) - |001) + |000), 

onto the orthogonal complement X"*", span a complex vector space of dimension six. ■ 



The following is now a direct consequence of Proposition [2] and Fact (TO) 



Theorem 2. The momentum map fiyy separates all K -orbits inside Xw. That is, 
every quantum state in the W SLOCC class is (up to local unitaries) uniquely determined 
by the spectra of its reduced density matrices. 



Points in the interior of ^{X\y) are also in the interior of ^(P('H)) and hence for any point 
a inside the interior ^{Xw) C ^(P(?^)) there is a point x G Mmax such that a = fi{x). 
The manifold n Mmax is 5-dimensional. Moreover, fi~^{fi{^a)) Pi M^^x is an 11- 

dimensional manifold consisting of 9-dimensional orbits K.y with the property fj,(K.y) = 
Qa- Every manifold H Mmax contains a unique orbit K.x with x £ G.x\y. We 

have illustrated the situation in Figure [2j 
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(c) The boundary of the Kirwan polytope 

In our analysis at the beginning of this section we did only consider quantum states 
that are mapped into the interior of the Kirwan polytope. We will now prove that any pure 
quantum state of three qubits that is mapped to the boundary of the Kirwan polytope 
^{¥{T-L)) is, up to local unitaries, uniquely determined by the spectra of its one-body 
reduced density matrices. 

Let us therefore consider a quantum state x G P(^) that is mapped to the boundary 
of the Kirwan polytope and write 

fi{x) = (diag(-iAi, iAi), diag(-iA2, iA2), diag(-f A3, zAs)), (4.5) 

according to ()4.2p . Then the minimal eigenvalue pj of the reduced density matrix pj is 
given by Pj = ^ — Aj. We distinguish two cases: 



Non-degenerate case We shall first treat the case where the reduced density matrices 
of X all have non-degenerate eigenvalue spectrum; that is, fj.{x) is contained in the interior 
of the positive Weyl chamber. We may assume by symmetry that /-i(x) is contained in 
the face corresponding to the equation pi =P2 +P3 (cf. Fact[9|), i.e. 

-Ai + A2 + A3 = ^. (4.6) 

That is, ijl{x) is orthogonal to (annihilated by) the Lie algebra element iS, G t, where 

e = \ (diag(l, -1), diag(l, -1), diag(l, -1)) , 

and Fact [6] implies that is contained in the Lie algebra of the isotropy subgroup of x. It 
follows that X = [v\ for some eigenvector v GTi of with eigenvalue 1/2. By diagonalizing 
the action of ^ on ^ we find that 

x = [ci|000) + C2IIIO) +C3IIOI)] 



Pure states of three qubits and their single-particle reduced density matrices 



15 



for some constants Cj. Any such state is contained in Xw, the closure of the SLOCC 
class of the VF-state, as can be seen by applying |0) -H- 11) to the first subsystem and 
comparing with ()4.4p . We can therefore use Theorem [2] to conclude that x is determined 
up to local unitaries by the eigenvalues of its reduced density matrices. 

Degenerate case We will now treat the case where at least one of the reduced density 
matrices is maximally mixed. Without loss of generality, we may assume that /i(x) is 
contained in the face of the Kirwan polytope defined by pi = ^, i.e. Ai = 0. Note that we 
cannot apply the same reasoning as above, since the faces Xj = arise from intersecting 
fi(P{T-L)) with the positive Weyl chamber and not from the geometry of the momentum 
map. By the Schmidt decomposition, and up a local unitary on the first subsystem, 

x=[^{\o)'^\^) + \i)0\m 

for orthogonal vectors {(f>\tp) = 0. The two-qubit reduced density matrix P23 of any such 
state is a normalized projector onto the two-dimensional subspace of Cj? with basis 
vectors and \tp); conversely, any other choice of orthonormal basis gives rise to a local 
unitarily equivalent state. 

Let us denote by G(2, 4) the Grassmannian consisting of two-dimensional subspaces 
/C C (g) C^. Similarly to the case of the projective space, we can consider the tensor 
product action of G' = SL{2) x SL{2) and of its maximal compact subgroup K' = 
SU{2) X SU{2). A momentum map for the if'-action is given by 

fi' ■.G{2,4) 3 IC^ i (^P2 - ^I, P3 - ^I^ ee , t'= su(2) © su(2) , 

where p2 and denote the reduced density matrices of P23 = ^Pk.i the normalized 
projector onto the subspace /C. In view of Fact [TOl it suffices the establish sphericality of 
this Grassmannian with respect to the action of G': 

Proposition 3. The G' -variety G(2,4) is spherical. 

Proof. Let B' denote the Borel subgroup consisting of lower triangular matrices in 
G' . As in the proof of Proposition [2l we will show that there exists a point /C G G(2,4) 
at which dime 7a:(-S'./C) = dime G(2, 4) = 4 (noting that G(2,4) is smooth). It will be 
convenient to work with coordinates. Let us therefore consider the Pliicker embedding 

2 

G(2,4)^P(/\C^), Spanc{|0),|^)}^[|./.) A|^)]. 

The image of the subspace /C = Spanj^{| + +), | )}, where |ib) = "^(|0) ± |1)), is 

x= [| ++) A I - -)] = [|00) A |01) + |00) A |10) - |01) A |11) - |10) A |11)]. 
It follows that the tangent space T]q[B' .fC) is spanned by the vectors 

(£;2i (g)/)xoc |00) A 111) - |01) A |10), 
[I ®E2i)x(x\m) A 111) + |01) A |10), 
((^11 - E22) (S> /)xoc |00) A |01) + |10) A 111), 
(/© (^11 - E22))x oc |00) A |10) + |01) A 111), 
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which are orthogonal to x and hnearly independent over C. We conclude that the tangent 
space is of complex dimension four. ■ 

In summary, we have proved the following result: 

Theorem 3. Let x be a pure quantum state of three qubits such that ^{x) is contained 
in the boundary of the Kirwan polytope ^{¥{T-L)) (i.e. its eigenvalues satisfy at least one 
of the inequalities in Fact [P| with equality). Then x is up to local unitaries uniquely 
determined by i.e. by the spectra of its one-body reduced density matrices. 



5. Summary 

In order to determine which pure states of three qubits are up to local unitaries uniquely 
determined by the spectra of their reduced density matrices, we have analyzed the change 
of the structure of the fiber ^~^(a) as a varies in the Kirwan polytope \I'(P('H)). We 
have shown that M„ = '^~^{a)/ K is generically a two-dimensional space. For the points 
in the boundary of the Kirwan polytope, the situation is rather different: We have proved 
that in this case ^~^{a) / K is a single point, i.e. the dimension drops by two compared 
with the points inside the interior. We have therefore identified all JC-orbits that are 
uniquely determined by the spectra of the reduced density matrices. In addition, we have 
examined the Kirwan polytopes ^{G.Xj) for all six three-qubit SLOCC classes (i.e. for 
the closures of the orbits of the complexified group G = K'^) and their mutual relation. In 
particular, we have proved that states from the so-called W SLOCC class are up to local 
unitaries separated by ^, i.e. each ET-orbit inside G.xy/ is characterized by the collection 
of spectra of the one-qubit reduced density matrices. 

Interestingly, the drop of the dimension of on the boundary of the Kirwan 
polytope has the following counterpart in (Christandl et al. 2012): The probability density 
/(q) of the eigenvalue distribution of the reduced density matrices of a randomly chosen 
pure state of three qubits vanishes precisely on the boundary of the Kirwan polytope. 
Since it is well-known that /(a)=volMa for regular points of the momentum map 
(Duistermaat & Heckman 1982, Meinrenken & Woodward 1999), it is reasonable to 
wonder whether an analogous statement could hold more generally (the main result of 
(Meinrenken & Sjamaar 1999) suggests that this might in fact be true). 

We also noticed that the polytope '^{X^y) associated with the W SLOCC class is 
of the same dimension as the polytope ^'(Xghz) corresponding to the GHZ SLOCC 
class, whereas for bi-separable and separable states the Kirwan polytope is of strictly 
lower dimension. On the other hand, it is known that for three qubits only the W and 
GHZ SLOCC classes represent genuinely entangled states. This intriguing relationship 
suggests that by looking at the polytopes corresponding to different SLOCC classes one 
can decide to what extent states from this class are entangled. We believe that this is 
not a coincidence and conjecture that similar phenomena should be present for A^-qubit 
systems, where > 3. 
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